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Abstract. Recent advances in interpreting the most accurate to-date measurement of atomic parity viola- 
tion in Cs are reviewed. The inferred nuclear weak charge, Qw( 133 Cs) = — 72.65(28) C x P t(36)thcor, agrees 
with the prediction of the standard model at 1 a level. Further improved interpretation is limited by an 
accuracy of solving basic correlation problem of atomic structure. We report on our progress in solving 
this problem within the relativistic coupled-cluster formalism. We include single-, double- and triple- elec- 
tronic excitations in the coupled-cluster expansion. Numerical results for energies, electric-dipole matrix 
elements, and hyperfine-structure constants of Cs are presented. 
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1 Introduction 

This year, year of 2006, marks 50 years of the fall of par- 
ity conservation, one of the hallmark discoveries of the 
20th century physics. The field has started with the sem- 
inal Lee and Yang paper [1] and discovery of parity viola- 
tion in the nuclear /3-decay [2]. Soon after this discovery, 
Zel'dovich [3] contemplated possibility of observing parity- 
nonconserving (PNC) signal in atoms. He concluded that 
the effect was too small to be of experimental significance. 
In 1974, however, M.-A. Bouchiat and C. Bouchiat [4] re- 
alized that the PNC is amplified in heavy atoms. They 
showed that the relevant PNC amplitude scales steeply 
with the nuclear charge Z. In atomic physics, the first 
P-violating signal has been observed in 1978 by Barkov 
and Zolotorev [5] in Bi atom. Over the following decades 
the experiments were refined, with PNC signal observed 
in several atoms. So far the most accurate measurement 
has been carried out in 133 Cs by Boulder group [6]. 

Rich history of atomic parity violation is examined in 
a number of review articles, e.g., Refs. [7,8] and a book [9]. 
Here, due to the space limitation, we restrict our attention 
to recent developments in interpreting P-violation in Cs 
atom and report on the progress towards more accurate 
interpretation of atomic PNC in this atom. 

The PNC measurements are interpreted in terms of the 
weak nuclear charge Qw, which quantifies the strength 
of the electroweak coupling between atomic electrons and 
quarks in the nucleus. The relation between Qw and the 



PNC amplitude, -Epnc, can be represented as 
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(1) 



where k is an atomic-structure factor. Apparently, the 
interpretation requires atomic-structure calculations of k 
with an accuracy that matches the experimental uncer- 
tainty in £pNO in particular, in Cs, the present theo- 
retical uncertainty dominates over experimental error bar 
resulting in an overall 0.6% uncertainty in the inferred 
value of Q w ( 133 Cs). 

The underlying theory of the electroweak interactions 
is well established and its predictions have been verified 
in a variety of experiments. Much of the present efforts 
are driven by searches for "new physics" beyond the stan- 
dard model (SM). We may distinguish between two ap- 
proaches to such searches: low- and high-energy experi- 
ments. Atomic parity violation probes the low-energy elec- 
troweak sector of the SM. While the relevant momentum 
transfer is just ~ 30 MeV, the exquisite accuracy of the 
interpretation constrains new physics at much higher ener- 
gies. Indeed, following arguments of Ref. [10], we find that 
the present 0.6% determination of the Qw( 133 Cs) probes 
the new physics at a mass scale of 
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where Gf is the Fermi constant. Apparently by reducing 
the uncertainty in the determination of the weak charge, 
|#Qvp|i we may access even higher mass scales. Such a re- 
duction in the uncertainty of Qw is the goal of our present 
efforts outlined in Sec. 3. 
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While discussing the impact of atomic PNC on parti- 
cle physics, it is worth noting that colliders are blind to 
certain "new physics" scenarios [10]. For example, atomic 
PNC is uniquely sensitive to extra Z bosons predicted 
in grand unified theories, technicolor models, SUSY, and 
string theories. Generally, atomic PNC is one of a few 
probes of electroweak coupling below the Z-pole. Com- 
peting and at the same time complementary experiments 
include determination of the weak charge of electron at 
SLAC [11] and ^-nucleus deep inelastic scattering by the 
NuTeV collaboration [12]. With an upgrade at Jefferson 
Lab, an even more accurate constraint on electroweak cou- 
pling is expected to come from a P-violating e — p scat- 
tering experiment (Qweak collaboration) [13]. Still, the 
atomic PNC will continue to play an important comple- 
mentary role, as the atomic PNC is sensitive to electron- 
neutron couplings, while the Qweak experiment will probe 
electron-proton coupling, i.e., a different linear combina- 
tion in terms of the up and down quark-electron couplings. 



2 Deviation from the Standard Model and 
recent theoretical progress 

Parity violation in Cs has been a subject of a recent con- 
troversy. To infer the value of the weak charge, one needs 
to combine measurements with atomic-structure calcula- 
tions. The uncertainty in the value of the inferred Qw is 
usually determined by summing experimental and theo- 
retical uncertainties in quadrature, 



a = {^oxpt + Otheor} 



1/2 



(2) 



In 1997, the most accurate to date (0.35%) measurement 
of PNC has been carried out by Boulder group [6] . At that 
time, the accuracy of calculations [14, 15] has been esti- 
mated to be 1%. In 1999, Bennet and Wicman [16], in light 
of new lifetime data which improved theory-experiment 
agreement, have reduced the theoretical error bar to 0.4%. 
While compared with the prediction of the SM, the result- 
ing value of Qw differed by 2.5 a, one of the largest devia- 
tions in the recent history. The deviation has prompted a 
substantial interest from particle physics community (see, 
e.g., [10, 17, 18, 19]). At the same time the reduced theoret- 
ical uncertainty raised the questions whether some " small" 
sub-1% atomic-structure effects could be the reason for the 
deviation. 

Over the last few years, there has been an impor- 
tant progress in understanding "small" corrections, such 
as Breit, radiative, and neutron-skin corrections. These 
advances are summarized in Table 1. 

Taking into account the tabulated "small" corrections, 
revised value of vector transition polarizability and some- 
what improved value of -Epnc [31], we arrive at the fol- 
lowing revised value of the weak charge (the assigned the- 
oretical uncertainty of 0.5% is consistent with an estimate 
of Ref. [31] and the central value with that of Ref. [29]) 



This value agrees with the SM value [19] of -73.19(13). 
at 1 a level. 

The present theoretical uncertainty in PNC calcula- 
tions is estimated to be 0.5%. Since the overall error is 
obtained with quadrature (2), reducing crthcor below the 
experimental uncertainty of 0.35% will either reveal new 
physics or set important constraints on competing exten- 
sions to the SM. In addition, it will set a stage for further 
experimental efforts, such as those by Paris group [32,33]. 
Other experimental efforts which will benefit from the im- 
proved calculations are PNC measurements with Ba + [34, 
35], Tl [36], and Fr [37]. 



3 High-accuracy atomic-structure calculations 

The present accuracy in the determination of the weak 
charge of Cs is limited by the accuracy of solving the ba- 
sic correlation problem. The many-body problem is hard. 
Even classically, the three-body problem can not be solved 
in closed form. While adequate numerical approaches (e.g., 
configuration interaction and variational methods) were 
developed for few-electron atoms, accurate solution of the 
many-body problem for atoms with many electrons still 
remains a challenge. 

In a particular case of the PNC problem in Cs, we need 
to evaluate the PNC amplitude for the 651/2 - » 7 £1/2 
transition 



Epnc = ^2 

+ E' 



(7S|J>|nP 1/2 )(nP 1/2 |g w |6S) 
E e s - E nPl/2 
(7S\H w \nP 1/2 )(nP 1/2 \D\6S) 



E7.S — E nPl 



(4) 



/2 



Here D and Hy/ are electric-dipole and weak interaction 
matrix elements, and Ei are atomic energy levels. The 
effective weak interaction, averaged over quarks is accu- 
mulated inside the nucleus 



w 



Gf 



Qw 75 p{r) 



(5) 



QwC^Ca) = -72.65(28) oxpt (36) theor 



(3) 



where p(r) is the neutron-density distribution. One could 
easily recast Eqs. (4) and (5) in terms of the structure 
factor k of Eq. (1). 

We wish to evaluate accurately the sum (4). Let us 
comment on pre-requisites for such calculations, (i) Due to 
the particular characters of the involved operators in sum- 
mation (4), the approximate wave- functions have to repro- 
duce the short-range (close to the nucleus) and long-range 
atomic properties simultaneously. The accurate short-range 
description requires ab initio relativistic approach, as the 
ratio of velocity of electron to the speed of light near the 
nucleus is ~ aZ <~ 0.4 for Cs. (ii) With respect to the ac- 
curacy, we notice that simple Dirac-Hartree-Fock (DHF) 
calculations for the hyperfine-structure (HFS) constant of 
the ground state are about 50% off from the experimental 
value. This constant describes strength of coupling of the 
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electron to the nuclear magnetic moment and its short- 
distance nature mimics behavior of the weak matrix ele- 
ments. Therefore we arrive at a conclusion that so-called 
correlation corrections (those beyond the DHF approxi- 
mation) have to be addressed. We treat the correlations 
within the systematic and successful methods of many- 
body perturbation theory. 

How do wc estimate the accuracy of calculating -Epnc? 
Ultra-precise experimental data are available for Cs. These 
data have accuracy better than our anticipated theoret- 
ical uncertainty thus allowing us to calibrate our calcu- 
lations. Atomic energies are known to many significant 
figures. Ultra-precise (0.01% accurate) value of lifetime of 
the 6P 3 / 2 state of Cs has been recently inferred from pho- 
toassociation spectroscopy [38,39] with ultracold atoms. 
The HFS constants for Cs are also measured to a good 
accuracy (the HFS constant A$s of 133 Cs is known ex- 
actly by the definition of the second) . 

To reiterate, to further improve interpretation of atomic 
parity violating signals, we need a high-accuracy ab initio 
relativistic many-body method capable of reaching the 
accuracy level of 0.1% for Cs. 



4 Relativistic coupled-cluster method 

Many-body perturbation theory (MBPT) provides a sys- 
tematic way of treating correlation corrections, the critical 
issue for an improved interpretation of atomic PNC. First 
we review main ideas of MBPT and the all-order coupled- 
cluster (CC) method. Then we describe our present CC- 
inspired computational scheme and illustrate it with nu- 
merical results for basic atomic properties of Cs atom. 



4.1 Generalities 

In MBPT the atomic Hamiltonian is partitioned as 

where ft. nuc includes the kinetic energy of an electron and 
its interaction with the nucleus, Udhf is the DHF poten- 
tial, and the last term represents the residual Coulomb in- 
teraction between electrons. In MBPT the first part of the 
Hamiltonian is treated as the lowest-order Hamiltonian Hq 
and the residual Coulomb interaction as a perturbation. 
The perturbative expansion is built in powers of residual 
interaction and the derivations typically involve methods 
of second quantization and diagrammatic techniques. 

One of the mainstays of practical applications of MBPT 
is an assumption of convergence of series in powers of the 
residual interaction. Sometimes the convergence is poor 



and then one sums certain classes of diagrams to "all or- 
ders" using iterative techniques. One of the most pop- 
ular all-order methods is the coupled-cluster (CC) for- 
malism [40,41]. It is widely employed in atomic and nu- 
clear physics, and quantum chemistry [42]. The relativistic 
atomic-structure CC-type calculations were carried out, 
for example, in Refs. [43,44,45,46,47,48]. 

The Hamiltonian in the second quantization (based on 
the DHF basis) reads 

H = H a + G 

= J2 £ * { a l a *} + 2 ^Isijkl \a\a)aia k } , (6) 

i ijkl 

where Hq is the lowest-order Hamiltonian and the resid- 
ual Coulomb interaction G is treated as a perturbation. 
The operators ai and a\ are annihilation and creation op- 
erators, and {...} stands for a normal product of operators 
with respect to the core quasivacuum state |0 C ). In the low- 
est order the atomic wave function with the valence elec- 
tron in an orbital v reads l^P^} — aj|0 c ). For example, v 
can represent the ground state 6s orbital of Cs atom. For- 
mally we can introduce a wave operator fl that produces 

the exact many-body wave function, \%) = Q\^ 0) ). 

Central to the CC method is realization that the wave 
operator can be represented via exponential ansatz [49] 

Q = { C MK)} = 1 + K + ± {K 2 } + ... . (7) 

The operator K can be compellingly separated into clus- 
ter operators combining simultaneous excitations of core 
and valence electrons from the reference state |^ r i°' ) ) to all 
orders of MBPT, 

K = S + D + T+--- , (8) 

i.e., K is separated into single (S), double (D), triple (T), 
and higher-rank excitations. For the univalent systems it 
is convenient to subdivide cluster operators into core and 
valence classes 

K = K C + K V . (9) 

Clusters K c involve excitations from the core orbitals only, 
while K v describe simultaneous excitations of the core and 
valence electrons. Then S — S c + S v , D = D c + D v , 
etc. The cluster operators satisfy so-called Bloch equa- 
tions [49]. 

While the CC formulation is exact, in practice the full 
cluster operator K is truncated at a certain level of exci- 
tations. If we restrict it to single and double excitations 

K = K^) ~ Sc + D c + S v +D v (10) 

we arrive at the widely employed coupled-cluster single 
double (CCSD) method. 

Linearized version of the CCSD method was employed 
by the Notre Dame group for high-accuracy calculations 
of various atomic properties [43,44,47,48]. In this approx- 
imation, one discards nonlinear terms in the expansion 
of exponent in Eq. (7), Q w 1 + if (SD) . We will refer to 
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this approximation as the singles-doubles (SD) method. 
The resulting SD equations are written out in Ref. [43] . A 
typical ab initio accuracy attained for properties of heavy 
alkali- metal atoms is at the level of 1%. 

Since the present goal is to reduce theoretical uncer- 
tainties to the level of 0.1-0.2% we have to go beyond 
the SD approach. A systematic step in improving the SD 
method would be an additional inclusion of triple and non- 
linear double excitations. However, considering the present 
state of available computational power, the full incorpo- 
ration of triples (specifically, core triples) seems to be un- 
managable for heavy atoms. For instance, for Cs storing 
and manipulating core triple amplitudes would require ~ 
100 Gb of memory. 

To motivate next-generation formalism, we have ex- 
plicitly computed 1648 fourth-order diagrams for matrix 
elements that appear due to triple excitations and non- 
linear terms (i.e., those omitted in the SD method) [50, 
51]. We observe from numerical results for electric-dipolc 
matrix elements in Na [51] and Cs [52] that the contribu- 
tions from valence triples T v and nonlinear doubles D n i 
are much larger than those from core triples T c . This 
leads to our present level of approximation: we discard 
core triples and core nonlinear terms and incorporate the 
valence triples and valence nonlinear terms into the SD 
formalism. The resulting approximation will be referred 
to as CCSDvT method. 



4.2 Driving equations in the CCSDvT approximation 

Below we write down the CC equations for cluster am- 
plitudes in the CCSDvT approximation. Here we present 
topological structure of the equations only. A detailed tab- 
ulation of the formulas can be found in our paper [53]. 
The equations in the SD approximation are presented in 
explicit form in Ref. [43]. The CCSDvT equations for the 
core cluster amplitudes S c and D c are the same as in the 
SD approximation. 

For valence triple amplitudes we obtain symbolically 

-[H ,T V ] + SE V T V « T V [D C ] + T V [D V ] . (11) 

Here [H ,T V ] is a commutator, and 5E V is the correlation 
valence energy defined as 

6E V = SE SD + 5E CC + 8E vT , (12) 

where correction SEsd is obtained within the SD approach, 
correction SEcc comes from nonlinear CC contributions 
and 5E v t is due to valence triples. Contributions T V [D C ] 
and T V [D V ] denote the effect of core and valence doubles 
on valence triples, respectively. At present we include only 
these effects omitting the effect of valence and core triples 
on valence triples (T„[T„] and T W [T C ]) and nonlinear CC 
contributions. These arc higher-order effects which com- 
putationally are much more demanding. 



The topological structure of the valence singles equa- 
tion is 

-[H ,S V ]+SE V S V w SD+ 
S V [S C ® S v ] + S V [S C ®S C ] + 
S v [S c ® D v ] + S v [S v ® D c ] + S v [T v ] . (13) 

Here S V [S C ® S v ] stands for a contribution from the ex- 
citations of core and valence electrons resulting from a 
product of clusters S c and S v . All other terms are defined 
in a similar fashion. 

Finally, equation for valence doubles can be symboli- 
cally represented as 

-[H ,D V ]+SE V D V w SD+ 
D V [S C <g> S v ] + D V [S C ® S c ] + 

D v [S c ®D V ]+D V [S v ® D c ] + D v [S c ® D c ] + (14) 
D v [D c ® D v ] + D v [S c ® T v ] + D v [S v ® T c ] + D v [T v ] . 

Solution of the above equations provides us with the 
cluster amplitudes and correlation energies. Numerical re- 
sults for the energies will be presented in Section 4.4. At 
this point, with the obtained wave functions we proceed 
to evaluating matrix elements. 

4.3 Matrix elements 

The SD method has already proven to be successful in 
calculations of various atomic properties. For heavy alkali- 
metal atoms the attained level of agreement with experi- 
mental data for the hyperfine constants is at 5% and the 
accuracy of a similar calculation for the electric-dipole am- 
plitudes is 0.5% (see, e.g., [48]). At the same time the 
accuracy required for our goals should be at the level of 
0.1-0.2%. In order to improve the overall accuracy we de- 
velop the technique of relativistic calculations of matrix 
elements beyond the SD approach. 

Given two computed CCSDvT wave functions, we may 
evaluate matrix elements of one-electron operator Z as 

wv ~ rrr — r~F rr? rr~\ ■ \ ) 

V(&w\&w)(1 r v\1 r v) 

The explicit expressions are given in Ref. [53]. Compared 
to the SD approximation, we include contribution of va- 
lence triples T v . They contribute both directly via ex- 
plicit contributions to matrix element formula and indi- 
rectly through modification of the SD amplitudes. 

It is worth pointing out, that the importance of the va- 
lence triples has been realized earlier by the Notre Dame 
group [44] . They have shown that at the SD level, the error 
for the HFS constants is as large as 5% for Cs. To rectify 
this problem, they proposed and implemented a scheme 
that approximates the effect S V [T V ], i.e., effect of valence 
triples on valence singles. While improving agreement for 
the HFS constants, their approximation leads to poorer 
agreement (compared to SD method) for the dipole ma- 
trix elements. The advantage of the Notre Dame scheme is 
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that it avoided expensive storing of triple excitations. Due 
to improved computational resources, we are able to store 
triples. Accounting for the triples in a rigorous fashion 
leads to a better agreement between theory and experi- 
ment. 

Compared to the Notre Dame approximation we also 
include dressing of matrix elements based on the CC ansatz. 
The idea of our method [5 2] is as follows. When the CC ex- 
ponent is expanded in Eq. (15), we encounter infinite num- 
ber of terms. We devised a method of partial summation 
(dressing) of the resulting series. Our formalism is built 
upon an expansion of the product of cluster amplitudes 
into a sum of n-body insertions. We considered two types 
of insertions: particle (hole) line insertion (line "dressing") 
and two-particle (two-hole) random-phase-approximation- 
like insertion. We demonstrated how to "dress" these in- 
sertions and formulated iterative equations. 

Another formal improvement over Notre Dame calcu- 
lations comes from including the CC nonlinear terms in 
the equations for valence singles (13) and doubles (14). 
We also include contribution of the core triples to matrix 
elements from direct fourth-order calculation. 

4.4 Numerical results 

Our developed numerical CCSDvT code is an extension 
of the relativistic SD code [47] which employs B-spline 
basis set. This basis numerically approximates complete 
set of single-particle atomic states. Here we use 35 out 
of 40 positive-energy basis functions. Basis functions with 
'max < 5 are used for singles and doubles. For triples we 
employ a more limited set of basis functions with l max (T v ) < 
4. Excitations from core sub-shells [4s,...,5p] are included 
in the calculations of triples while excitations from sub- 
shells [Is,.. .,3d] are discarded. 

Table 2. Contributions to removal energies of 6s, 6pi/ 2 , 
and 6P3/2 states for Cs in cm -1 in various approximations, 
^extrapolated correction is obtained by computing SD proper- 
ties with increasingly larger basis sets and interpolating them 
to I = co [54]. A comparison with experimental values is pre- 
sented in the lower panel. 





6s 


6pi/2 


6f>3/2 


-Bdhf 


27954 


18790 


18389 


5Esr> 


3869 


1611 


1623 


<5i?CCSDvT 


3350 


1387 


1220 


rptot 

^CCSDvT 


31304 


20178 


19608 


QED 1 


18 


-0.4 





A TP 101 

u ^extrapolated 


30 


20 


20 


-^filial 


31352 


20198 


19628 


-S'cxperim 


31407 


20228 


19675 



Reference [28]; Reference [55]. 

Computed removal energies of 6s, 6pi/ 2 , and 6p 3 / 2 
states of atomic cesium are presented in Table 2. The dom- 



inant contribution to the energies comes from the DHF 
values. The remaining (correlation) contribution is given 
by Eq. (12). Wc computed this correlation correction in 
SD and CCSDvT approximations. As it follows from the 
table the agreement with experiment is at the level of 0.1- 
0.2% for all considered states. We anticipate that including 
other corrections missed at this stage (e.g., nonlinear cor- 
rections to core amplitudes and core triples) can further 
improve the agreement with the experimental results. 

With the computed wave functions of the 6s, 6pi/2 and 
6P3/2 states we determine matrix elements. Numerical re- 
sults for magnctic-dipole hyperhnc-structure constants A 
and electric-dipole transition amplitudes are presented in 
Table 3. This Table is organized as follows. First we list 
the DHF and SD values, and the differences between CCS- 
DvT and SD values, Z\(CCSDvT) = CCSDvT - SD. We 
base our final ab initio results on the most sophisticated 
CCSDvT values. These values also include all-order dress- 
ing, and corrections due to core triples, computed in the 
the fourth order of MBPT. QED corrections are included 
where available. The results for the HFS constants include 
finite-nuclear size (Born-Weisskopf) effect. 

We find an excellent, 0.1%-level, agreement for dipolc 
matrix elements and the HFS constant of the ground state. 
The agreement for the the HFS constant of the 6p\/2 is 
only at 1% level. We are presently working on testing sen- 
sitivity of this constant to higher-order effects. 

We would like to emphasize that presently the correla- 
tion corrections at the level of a few 0.1% are comparable 
to radiative corrections. In this regard it would be use- 
ful to carefully compute them to unmask the remaining 
many-body effects. 



5 Summary and outlook 

Atomic parity violation plays an important role in testing 
low-energy electroweak sector of standard model. Interpre- 
tation of experiments in terms of nuclear weak charge re- 
quires calculations challenging the capabilities of modern 
atomic theory. Over the last few years, we have witnessed 
a substantial advance in evaluating corrections to parity 
violating amplitudes in heavy atoms. These small (sub- 
1%), but important corrections include Breit, radiative 
(vacuum polarization, self-energy, and vertex) and neu- 
tron skin corrections. As a result of this progress, the most 
accurate to date measurement of atomic parity violation 
in Cs has been brought into substantial agreement with 
the prediction of the standard model. 

Presently the theoretical interpretation is clouded by 
uncertainties in solving the basic correlation problem of 
atomic structure. In this paper we outlined our next-gen- 
eration many-body formalism for solving this problem. 
We tested our coupled-cluster-inspired method by com- 
puting basic atomic properties of Cs atom. All the com- 
puted properties are important for quantifying accuracy of 
the calculations of parity-violating amplitudes. Wc find an 
agreement at 0.1% for the ground state hyperfine struc- 
ture constant, El transition amplitudes, and energies. 
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However, a relatively poor 1% agreement of the HFS con- 
stant A for the 6P1/2 state with experiment requires fur- 
ther improvements of the method. The advantage of the 
employed coupled-cluster method is that it allows for such 
systematic improvements. It is anticipated that the further 
theoretical progress will refine constraints on new physics 
beyond the standard model and enable next round of ex- 
perimental studies. 

This work was supported in part by the US National Science 
Foundation, by the US NIST precision measurement grant, and 
by the Russian Foundation for Basic Research under Grant 
Nos. 04-02-16345-a and 05-02-16914-a. This manuscript was 
completed with a partial support by NSF through a grant for 
the Institute for Theoretical Atomic, Molecular, and Optical 
Physics at Harvard University and Smithsonian Astrophysical 
Observatory. 
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Deviation (1999) 
Bennet & Wieman [16] 


2.5a 


Based on calculations by 
Dzuba et al. (1989) [14] and 
Blundell et al. (1990) [15]. 


Theoretical correction 


Pull 


Reference 


Breit interaction 


— 1.2cr 


Derevianko (2000) [20]; 
Dznha pi al C20011 1211- 
Kozlov et al. (2001) [22]. 


Vacuum polarization ( +0.8a) 
Vertex/self-energy ( —1.3a) 

Neutron skin 


-0.5(T 

-0.4a 


Johnson et al. (2002) [23]; 
Kuchiev & Flambaum (2002) [24]; 
Milstein et al. (2002) [25,26]; 
Sapirstein et al. (2003) [27]; 
Flambaum & Ginges (2005) [28]; 
Shabaev et al. (2005) [29]. 
Derevianko (2002) [30] 



Updated correlated Bpnc and 

e — e P- violating interaction, 
Renormalization Q — > 0, 
Giant nuclear resonance 



+0.7a Dzuba et al. (2002) [31]. 



-0.08a 



Sushkov & Flambaum (1978) 
Milstein et al. (2002) [25] 



Total deviation (agreement) 



la 



Table 1. Present status of the 2.5a deviation of inferred Qw ( 133 Cs) from the prediction of the Standard Model. In the bulk of 
the Table, we summarize the recent theoretical progress on computing "small" sub-1% corrections to the 6s — 7s PNC amplitude. 
These corrections directly affect the inferred value of the weak charge of Cs. For illustrative purposes all the corrections (pulls) 



are tabulated in terms of the original Bennet and Wieman [16] uncertainty, a = { a 2 



1/2 



0.5%. 



Table 3. Magnetic-dipole hyperfine structure constants A (in MHz) and matrix elements of electric dipole moment (in a.u.) for 
133 Cs. Results of calculations and comparison with experimental values are presented. See text for the explanation of entries. 





A(6s) 




(6pi /a ||£>||6«) 


<6P3/ 2 ||£>| 


DHF 


1425.4 


160.94 


5.2777 


7.4264 


SD 


2438.0 


310.71 


4.4829 


6.3075 


Zi(CCSDvT) 


-136.9 


-20.92 


0.0256 


0.0363 


Complementary 


corrections 






Line dressing 


-12.5 


-2.16 


0.0094 


0.0107 


Vertex dressing 


4.3 


0.29 


-0.0067 


-0.0088 


MBPT-IV (core triples,...) 


7.8 


1.14 


0.0001 


0.0001 


Breit + QED 


-6.5 1 




0.0024 2 




Extrapol. for I — oo 3 


5.0 


0.37 


-0.004 


-0.006 


Final CCSDvT + corrections 


2299.2 


289.43 


4.5097 


6.3398 


Experiment 


2298.2 


291.89(8) 4 


4.5049(17) 5 


6.3404(3)' 



^ef. [56]; 2 Refs. [20,28]; 3 Ref. [54]; 4 Ref. [57]; 5 Refs. [58,39]; 6 Ref. [39]. 



